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ABSTRACT 

We investigate scaling relations between the dark matter (DM) halo model pa- 
rameters for a sample of intermediate redshift early - type galaxies (ETGs) resorting 
to a combined analysis of Einstein radii and aperture velocity dispersions. Modeling 
the dark halo with a Navarro - Frcnk - White profile and assuming a Salpcter initial 
mass function (IMF) to estimate stellar masses, we find that the column density S 
and the Newtonian acceleration within the halo characteristic radius r s and effective 
radius Reff are not universal quantities, but correlate with the luminosity Ly, the 
stellar mass M* and the halo mass A^ooi contrary to recent claims in the literature. 
Wc finally discuss a tight correlation among the DM mass Mum (Reff) within the 
effective radius Reff, the stellar mass M+(R e ff) and Reff itself. The slopes of the 
scaling relations discussed here strongly depend, however, on the DM halo model and 
the IMF adopted so that these ingredients have to be better constrained in order to 
draw definitive conclusions on the DM scaling relations for ETGs. 
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The current cosmological paradigm, the concordance ACDM 
model, re lies on two main co mponents, namely dark en- 
ergy (e.g. ICarroll et al. 1 1 19921 ) and dark matter. Although 
in excellent agreement with all the cosmological probes 
|Komatsu et al. 1120091 ; IPercival et al. 1120091 : lLampteil et al. I 
120091 ). the ACDM model is nevertheless afflicted by serious 
problems on galactic scales. In this framework, the forma- 
tion of virialized DM haloes from the initial tiny density 
perturbations is follo wed at later stages through numerical 
N-body simulations (jBertschinger Ill998l ). It became appar- 
ent that the spherically averaged density pro file, Pdm(t), of 
DM haloes is independent of the halo mass (jNavarro et al. I 
1 19971 ) and well described by a double power - law relation 



galaxies : kinematic and dynamics - galaxies : elliptical 



better constrain the formation scenario and the DM proper- 
ties. Recently, much work has b een dedicated t o this issue. 
Using a sample of local ETGs, iTortora et al~l (|2009l . T09 
hereafter) have found that DM is t he main driver of the 
Fundamental plane (FP) tilt (see also Cappellari et al.ll2 006. 



Bolton et al.ll200lmvde fc Bernardill2009l . lGraves fc Faber I 



2010) . lAuger et al. 1 12010a ) and that the average spherical 



with pdm oc 



in the outer regions and pdm oc 



with a > centrally. On the contrary, observations of spiral 
galaxies seem to de finitely point t owards cored models, i.e. 
a = at the centre (|de Blok 1120091 ). Understanding whether 
such a discrepancy is due to some physical process not cor- 
rectly modeled in simulations or to a failure of the CDM 
paradigm is still a hotly debated issue. As a valuable tool 
to address this problem, one can look for scaling relations 
among DM halo parameters and stellar quantities in order to 



DM density is a dec reasing function of stellar mass (see also 
iThomas et al.ll2009l ). Based on data from rotation curves of 
~ 1000 spiral galaxies, the mass models of individual dwarf 
and spiral g alaxies and the weak lensing sig nal of elliptical 
and s pirals, iDonato et al. I (|2009l . D09) and iGentile et al. I 
(J2009I . G09) have found strong evidence for the constancy of 
the central D M column density over 12 orders of magnitude 
in luminosity. iNapolitano. Romanowskv fc Tortora I (|2010l . 
NRT10) have shown that, on average, the projected den- 
sity of local ETGs within effective radius is systematically 
higher than the same quantity for spiral and dwarf galaxies, 
pointing to a systematic i ncreas e with halo mass as sug- 
gested bv lBovarskv et al. I (|2009al . B09), who have extended 
the samples analyzed above to both group and cluster scale 
systems. 

In order to try to discriminate between these contrast- 
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ing results, we present here an analysis of the DM scal- 
ing relations for a sample of ETGs at intermediate redshift 
({z) ~ 0.2) using lensing and velocity dispersion data to 
constrain their parameters. The mass models, the data and 
the fitting procedure are described in fj2] In ij3]we describe 
the main results, while (JUis devoted to a brief review of the 
results and conclusions. 



2 ESTIMATING MASS QUANTITIES 

As a preliminary mandatory step, we need to determine the 
quantities involved in the above scaling relations. To this 
end, one has first to choose a model for the stellar and DM 
components and then fit the observational data in order to 
infer the quantities of interest from the constrained model. 



2.1 Stellar and DM profiles 

Motivated by the well known result that the surface 
brightness profiles of ETGs are well fitted by the Sersic 
(1968) law, we des c ribe the stellar component with the 
Prugniel fc Simien I (119971 . PS hereafter) profile (see also 
Marquez et al. Il200ll for further details). The choice of the 
DM halo model is quite controversial. Rotation curves of 
z = spiral galaxies are better fitted by cored models 
(|de Blok II2009P ). but we are here considering ETGs at inter- 



mediate z so that it is not straightforward to extend these 
results to our case. In order to explore th e impact of the 
DM halo profile, we therefore adopt both a lBurkert I (| 19951 ) 
model with 



Pbt\ 



PB{r) (r + r B )(r 2 + r|) ' 

and an NFW (|Navarro et al. II1997T ) profile with 



PNFwij-) 



p s r s 



r(r + r s ) 2 



(1) 



(2) 



Both the 3 - dimensional and projected masses M(r) and 
Mp ro j (R ) can be analytically e valua ted, with Mp ro j(R ) 
given in iPark fc Ferguson I (|2003h and iBartehnannT (|l996h 
for the Burkert and NFW models, respectively. 

In order to constrain the model parameters, we rely on 
the estimate of the projected mass Me = M pro j(RE) in- 
ferred by the measurement of the Einstein radius Re in 
a lens system. While lensing probes the mass projected 
along the line of sight, t he aperture velocity dispersion a ap 
IjMamon fc Lokas 1 120051 s ) provides complementary informa- 
tion on the internal dynamics thus strengthening the con- 
straints. 



2.2 Data and fitting procedure 

We make use of the sample o f 85 lenses collected by the Sloan 
Lens ACS (SLACS) survey (|Auger et al. 1120091 ) and first se- 
lect only ETGs with available values of both the velocity 
dispersion a ap (measured within an aperture of i? ap = 1.5") 
and the Einstein radius Re, thus ending up with a dataset 
contai ning 59 objects. For each lens, we follow I Auger et al. I 
(|2009l ) setting the Sersic index n = 4 and the effective radius 
R e ff and total luminosity Lv to the values inferred from 
the V - band photometry. The SLACS collaboration has also 



provided an estimate of the total stel l ar ma ss (th eir Table 
4) fro m which we use both IChabrierl (j200lh and ISalpeterl 
(|l955h IMFs to investigate the effect of the IMF on the scal- 
ing relations. 

Before fitting the model to the data, it is worth noting 
that, for both the Burkert and NFW models, the halo pa- 
rameters are different from one lens to another. Since we 
have only two observed quantities, namely (Me, (Tap) for 
each lens, determining (px,rx) (with X = B or s for the 
Burkert and NFW models, respectively) on a case -by -case 
basis would give us very weak constraints. We therefore bin 
the galaxies in 10 equally populated luminosity bins (the last 
one actually containing 5 objects) and resort to a different 
parametrization using quantities that it can be reasonably 
assumed to be the same for all the lenses in the same birjj. 
As one of the fitting parameters, we choose the virial M/L 
ratio, T v i r = M V i r /Lv with M V i r the DM halo mass at the 
virial radiu^J R v i r . Should T V i r depend on Lv, our luminos- 
ity bins are quite narrow so that any change in T v i r should 
be so small that can be safely neglected. We then use log r\x 
with r\x = Tx /Reff as the second parameter. 

Note that, although (T„; r ,log?7x) are the same for all 
the galaxies in a bin, (M V i r ,rx) still change from one lens 
to another thus allowing us to estimate all the quantities 
of interest on a lens -by -lens basis. In order to constrain 
(T v i r , log rix), we maximize a suitable likelihood function, 
composed of two terms, the first (sec ond) one referring to 
the lensing (dynamics) constraints (see lCardone et al. 1120091 
for further details) . In order to efficiently explore the param- 
eter space, we use a Markov Chain Monte Carlo (MCMC) al- 
gorithm running chains with 100000 points reduced to more 
than 3000 after burn in cut and thinning. For a given galaxy, 
we compute the quantities of interest for each point in the 
chains and then use Bayesian statistics to infe r median val- 
ues an d 68% confidence intervals. We follow ID' Agostini I 
(J200J) to correct for the asymmetric errors in our estimates. 



3 RESULTS 

Before investigating scaling relations, it is worth checking 
whether our PS + DM model works in fitting the lens data. 



3.1 The fiducial stellar + DM profile 

We have considered four PS + DM models by combining 
the two halo profiles (Burkert or NFW) with the two IMFs 
(Chabrier or Salpeter) adopted to set the total stellar mass. 
We find that all the four models fit remarkably well the 
lensing and dynamical data with rms deviations ~ la and 
the lens observed values deviating no more than 2a from the 
model ones for most of the bins. On the one hand, such a re- 
sult ensures us that the estimates of the different quantities 
entering the scaling relations we are interested in are based 
on empirically motivated models. On the other hand, this 



1 Note that this procedure allows us to fit for these two alter- 
native parameters on a bin - by - bin basis thus having 2J\f\ >in — 2 
degrees of freedom with Afbin the number of lenses in the bin. 

2 We define R v i r as the radius where the mean density 
M v i r /(4:/3)nR 3 . equals A c (z)/jj.f(z) with A c (z) as in Bryan & 
Norman (1998) and pM the mean matter density at z. 
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Table 1. Constraints on halo parameters for the Burkcrt models with Chabrier and Salpeter IMF. The median luminosity from the 
galaxies in each bin is reported in column 1, while the adjacent columns give the maximum likelihood parameter (T„; r ,log?js)Ml,, the 
value of — 21n£ at the maximum, median value and 68 per cent confidence interval for the DM parameters (T v i r , log»7 s ) for each model. 
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Table 2. Same as Table 1, but for NFW models. 
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same result tells us that the data we are using are unable 
to discriminate among different choices. This is an outcome 
of the limited radial range the data probe. Indeed, given 
typical values of ReJR^jj and R ap /R e ff, both Me and a ap 
mainly probe the inner regions of the lenses so that they can 
provide only weak constraints on the behaviour of the mass 
profile in the outer DM halo dominated regions. 

Although statistically equivalent (as can be quantita- 
tively judged on the basis of the close values of £ mal in 
Tables 1 and 2) , the four models may be ranked by examin- 
ing the constraints on their parameters. Indeed, the values 
in Table 1 show that the PS + Burkert models fit the data 
with quite low values of T V ir thus leading to unexpectedly 
small virial masses. Roughly averaging the maximum like- 
lihood T V ir of the different bins, we get (T V i r ) — 22 ± 14 
(12 ±9) M/L v ,q using the Chabrier (Salpeter) IMF. On the 
contrary, when adopting the NFW profile, we find (Y„j r ) = 
590 ± 210 (360 ± 150) M/L V , Q for a Chabrier (Salpeter) 
IMF. Our results from NFW + Salpeter are qualitatively 
consistent with previous estimates of T^r ~ 200 M/Ly.Q 
obtained relying on galaxy - galaxy lensing (jGuzik fc Seliak I 
120021 ; iHoekstra et ai. 1120051 ) o r combining strong a nd weak 
lensing and central dynamics (JGavazzi et al. 1120071 ). There- 



fore, we will consider the NFW + Salpeter as our fiducial 
modeo retaining the other cases just for investigating the 
dependence of the scaling relations on the halo profiles and 
IMF adopted. This assumption is also confirmed by recent 
findings which point to a Salpe ter IMF when a n uncon- 
tracted NFW profil e is assumed (JTreu et alj|2010l ; NRT10, 
lAuger et altolOal ). 

The constraints on (T„i r ,rog?7 s ) for the four cases con- 
sidered are summarized in Tables 1 and 2 for Burkert and 
NFW models respectively. The marginalized constraints on 
T v ir show that the NFW + Salpeter model provides reason- 
able values for the virial M/L ratio in all bins (but the first) 
thus motivating our choice as a fiducial case. The 68% CL 
for Tvir are, however, quite large and asymmetrically ex- 
tended towards very large values. This can be qualitatively 



3 In order to put this choice on firmer statistical grounds, we 
could have redefined our likelihood by adding a prior on Y„j r 
rather than the uninformative flat one we have used here. How- 
ever, since previous estimates of Y„i r are affected by large errors 
and based on model assumptions, we have preferred to avoid the 
risk of a fit driven by the prior examining a posteriori the resulting 
values. 
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Figure 1. S(r B ) and S(R e ff) for the NFW + Salpeter model as a function of stellar mass (left panels) a nd S( r s ) vs halo mass M200 
(right panel), shown as black points. We superimpose the best fit linear relation using the ID' Agostini I l|2005l ) (red solid) and direct 
fit (blue dashed) methods. On the right axes is shown the equivalent acceleration scale goM-, defined as in the text and derived by 
the fit of S — gDM relation. The black error bars set the median data uncert ainties. The res u lts w hen a Chabrier IMF is used are 
shown as open boxes. In the left -bottom panel, the median from NRT10 and iTortora et al. I 1120101 ) are shown as red symbols and 
gray region, respectively (in both the cases the la scatter of the distribution is plotted). In the right panel, we add the B09 best fit 
linear relation, the secondary - infall mo del (Boyarsky et al. 2009b) prediction and the results from the ACDM N - body simulation of 
iMaccio. Dutton fc van den Bosch] J200gt ). 



explained noting that the larger is T V i r , the larger is M v i r 
and hence R V i r . In such a case, Re/Rvit and R ap /Rvir be- 
come quite small such that the data are less and less able to 
constrain the outer regions. As a consequence, the marginal- 
ized likelihood function £(T„; r ) has a long flat tail for large 
T V i r thus giving rise to the reported asymmetric errors. A 
similar effect also explains why the errors on log rj 3 are asym- 
metric and still large. The wide confidence intervals prevent 
us from quantifying how Y„ir and log rj 3 change from one bin 
to another. Excluding the first bin, characterized by large 
uncertainties, we find that both {T V i r ,logrj s ) do not show 
any significant trend with Lv- 



3.2 DM correlations 

We now investigate the relation between the column den- 
sity S(R) — M pro j(< R)/irR 2 , luminosity, stellar and halo 
masqj for the lenses in our sample. The main results of this 
analysis are shown in Fig. [T] Despite the small mass range 
probed and the large errors, we find that S(r a ) is positively 
correlated with Lv, M* and M200, confirming the results in 
B0 9. The slope of the correlations may be estimated using 
the ID ' Agostini I (|2005l ) fitting method, which takes into ac- 
count the errors on both variables and the intrinsic scatter 
o~int- Concentrating on the trend with stellar mass, for the 
maximum likelihood parameter^, we get 



4 Note that, in order to be consistent with B09, we use A^oOi i-e. 
the mass within the radius -R200, where the mean density is 200 
times the cosmological mean matter density. 

5 We will refer always to the maximum likelihood values only, but 
the reader must be aware that, in a Bayesian framework, what is 
most important are the marginalized values and their confidence 



log S(r s ) = 0.29 log 



Af* 



10 11 M Q 



+ 2.17 



with an intrinsic scatter a% n t — 0.01. Although the results 
from a direct fid are only slightly different, for completeness 
we will plot the best fitted relations from both the methods 
(see Fig. [TJ. Similar results are found when the column den- 
sity is fitted as a function of Ly. 

If we fit the same relation, but replace S(r s ) with 
S(Reff), the trends are shallower, and the zeropoints change 
too. As shown in Fig. [l] an error weighted average over the 
galaxies in the sample gives {log S(R e ff)) ~ 3.1 in agree- 
ment, within the scatter, with the median logS{R e ff) ~ 2.9 
from the local ETC sample o f T09 and NRT10 and the 
results in ITortora et al. I (|2010r ) where a different analysis, 
based on an isothermal profile, is used on the same lens 
sample. All these results are in agreement with a scenario 
where ETGs surface from the merging of late - type systems 
so that, at a fixed radius, their DM content is larger than 
the one for spirals and dwarves (logS{R e ff) ~ 2 — 2.5, see 
NRT10 for further details). Should we use a Chabrier IMF 



ranges given in Table \3\ Since the median value of the slope and 
the scatter are close to the maximum likelihood ones, this choice 
has no impact on the discussion. 

6 The intrinsic scatter accounts for the deviations of the single 
galaxies from the underlying model leading to the fitted relation. 

7 In the direct fit, we minimize the usual \ 2 assuming that the 
errors on the x variable are negligible and no intrinsic scatter is 
present. These simplifying assumptions do not hold for most of 
the scaling r e lation s we have considered so that we resort to the 
ID' Agostini I (120050 method as a more reliable procedure. How- 
ever, if (Ti n t is small and the error on x lower than that on y, the 
two methods converge to the same maximum likelihood values. 
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Table 3. Marginalized constraints on the scaling relation parame- 
ters for the correlations involving the column density cS assuming 
the fiducial NFW + Salpetcr model. Columns are as follows: 1. 
correlation id; 2., 3., 4. median value and 68 confidence ranges 
for (or, 0, (Ti n t), where the linear relation logy = a logic + ft is 
fitted and (Tj„t is the intrisinc scatter. Note that, due to the fit- 
ting method, for each value of (a, Oj n j), is set by the condition 
that the likelihood is maximized, in other words we analytically 
marginalize over /3 when determining the maximum likelihood 
parameters. 

We warn the reader that, in the fit, we use the luminosity Ly 
in units of 10 11 L@ and the stellar (halo) mass in units of 
10 11 (10 12 ) Mq to reduce error covariance. 



id 



(«):£ 



wt\: 



[&int 



-Icr 



log S(r s )- log L v 0-28±g;i| 

log S(r 3 )- log Ah 0-29±g"_J| 

logS(r s )-logM 200 O.U+° f 

logS(U e// )-logM* O.Ut° \l 



2 34+0.01 
z -° -o.oi 

2 16+ - 10 
z - 1D -0.09 
n n+0.24 
z - x± -0.22 
2 OA+0-07 
z - 3 -0.08 



n o ? +0.043 
u.uo/_ 027 

o 7 +0.043 
u.uoi_ 027 

040+ ' 044 

n n57 +o.()39 
u.uo(_ 032 



(thus lowering the stellar masses by a factor ~ 1.8) we get a 
larger DM content and S(R e ff) which is a constant function 
of Mi, (see open boxes in Fig. [T]). 

A weaker positive correlation is found when plotting 
S(r 3 ) vs M200, the maximum likelihood fit being 



log 5(r s ) = 0.16 log 



M200 

10 12 Mr; 



2.11 



with tjint = 0.010. Our best fit relation is shallower than 
the B09 one, although the slope is consistent with their one 
(0.21) within the large error bars. Note that we have here 
explicitly taken into account the correlated errors on S(r s ) 
and M200 , while we do not know whether the fitting method 
adopted by B09 does the same. We therefore cannot exclude 
that the difference in slope is only an outcome of the use of 
different algorithms on noisy data. Actually, our S(r 3 ) val- 
ues are on average smaller than those in B09 (over the same 
mass range). Our results are also systematically smaller than 
the estimates fro m the ACDM N - body simulations of iso - 
lated haloes from iMaccio. Dutton fc van den Bosch l (|2008l ) 
and the prediction s from the secondary infall modelj (B09, 
iDel Popolo II2OO9I ). Although a wider mass range has to be 
probed to infer any definitive answer, we argue that the 
larger values in the literature are a consequence of neglect- 
ing the stellar component. The agreement with B09 would 
improve if a Chabrier IMF was used (see open boxes in the 
right panel of Fig. [TJ and, as discussed above, the trend 
would flatten out. 

When discussing the results for our reference model, 
the observed correlations argue against the universality of 
the column density proposed in DOQj. It is worth investigat- 
ing why we and D09 reach completely opposite conclusions. 



8 Note that the predictions for the secondary infall model of Del 
Popolo arc actually smaller than the B09b model ones (Del Popolo 
et al., in preparation). 

9 Actually, what D09 refers to as universal quantity is the prod- 
uct Pb^b- It is, however, possible to show that Snfw{ t s) — 
<Sg(1.6rg) (B09). It is then just a matter of algebra to get 
SbO-S^b) ~ 1.89pBra so that the constraint log(pBrs) = 



As a first issue, we note that D09 describe the DM halo 
adopting the Burkert model. Should we use this model to 
infer 5(1. 6vb), a Chabrier IMF (which does not strongly 
differ from the IMFs used in D09), and plot as a function 
of luminosity to be uniform with them, the best fit relation 
would have been 



log5(1.6r B ) = 0.02 log 



10 11 L 



2.65 



with aint = 0.15. In such a case, within a very good ap- 
proximation, we can assume 5(1. 6re) is indeed constant 
with Ly- An error weighted average of the sample values 
gives (log5(1.6r B )) = 3.07 ± 0.02 ± 0.20, where here and in 
following similar estimates, the first error is the statistical 
uncertainty and the second one the rms around the mean. 
Taken at face value, this estimate is significantly larger than 
the D09 one (similarly to what was found when discussing 
the comparison among the S(R e ff) values for ETGs and 
spirals), even if it is marginally within 2oy ms . However, as 
discussed above, the slope of the Sb{^-&tb)- Ly relation is 
strongly model dependent. Indeed, changing to a Salpeter 
IMF, we find 5B(1.6r_e) rx Ly 20 thus arguing in favour of 
non universality. Actually, the Burkert model is likely to 
be unable to fit the data for ETGs with reasonable values 
of Y v ir so that one should rely on the NFW + Salpeter re- 
sults to thus conclude that the column density 5 depends 
on mass for intermediate redshift ETGs. Investigating the 
reasons why NFW models are preferred over Burkert ones 
for ETGs, while the opposite is true for spirals can provide 
important constraints on galaxy formation scenarios, but it 
is outside our scope here. 

Motivated by the D09 findings on the universality of 
Pbtb and noting that the DM Newtonian acceleration, 
gDAi(r) — GMoM(r)/r 2 , is proportional to pbtb for a 
Burkert model, G09 have recast the D09 results in terms of a 
constancy of goMirs)- Then they extended this result show- 
ing that also the stellar acceleration, g*(r) — GMi,(r)/r 2 
at tb is a universal quantity. Using their same Burkert 
model and a Chabrier IMF, we indeed find that gDM(rs) 
is roughly constant with Lv, with an error weighted value 
(log </_da/(t"b)) = —8.89 ± 0.02 ± 0.20, in satisfactory agree- 
ment with the G09 result (loggrrjM (rs) — —8.5). Having 
said that, we find that g*(rs) is not constant, the best fit re- 
lation predicting <?*(t\b) oc L v ' 21 . We stress, however, that 
the situation is completely reversed if we use a Salpeter IMF 
giving gDM{rB) oc L v 20 and g*(rs) oa L y ° , so that draw- 
ing a definitive answer on which quantity is universal is an 
ambiguous task. Whatever is the correct IMF, we can never- 
theless safely conclude that the DM and stellar Newtonian 
accelerations cannot both be universal quantities, in con- 
trast with the claim in G09. Moreover, should the IMF be 
universal and intermediate between Chabrier and Salpeter, 
one could argue that neither gDM(rs) and g*(rs) are uni- 
versal quantities. 

As already discussed, our fiducial model is the 
NFW + Salpeter so that it is preferable to use this model 
when investigating whether the Newtonian accelerations are 
constant or not. Moreover, rather than using the halo char- 



2.17 ± 0.20 in D09 translates into log<S s (1.6r fl ) 
which we use as a a comparison value. 



2.44 ± 0.47 
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Table 4. Marginalized constraints on the scaling relation parame- 
ters for the correlations log gDM(Reff) -log Ly (upper part) and 
log g*(Reff) -log Ly (lower part) with the accelerations in m/s 2 
and the luminosity in units of 10 11 Lq. Columns are as follows : 1. 
model id, 2., 3., 4. median value and 68 per cent confidence ranges 
for (a, /3, <Ti n t). The model ids are BSC for Burkert + Chabrier, 
BSS for Burkert + Salpeter, NSC for NFW + Chabrier, NSS for 
NFW + Salpeter. Note that, since g*(R e ff) refers to stellar quan- 
tities only, the correlations involving this quantity are indepen- 
dent on the halo model. 



Model Id 



(« 



wtl: 



I \ + lo 

\<Tint)_ lo 



BSC 


-n i^+ - 145 

u.ioo_ 145 


-9 090+ 0003 

3.U3U_ 003 


_ +0.023 
u,1 '°-0.019 


BSS 


-0.004±S"iS 


Q 099+0. 004 
3.OZZ_ 002 


15fi+0° 33 
U.lOD_ 02g 


NSC 


0.088+ ;iii 


-9 926+ - 004 
y - aZO -0.004 


066+ 0043 
u.uou_ 035 


NSS 


0.245+°;^ 


QSA1 +0-007 


n 040+ - 043 

U.U4U„ 028 



NSC 
NSS 



-0.461 
-0.461 



+0.170 
-0.168 
+0.170 



-9 649+ - 004 

».Ot3_ 00 g 



-9.403 



+ 0.004 



0.178 
0.178 



+0.028 
-0.022 
+0.028 



acteristic radius (which refers to a different mass content 
depending on the model), we will discuss the results at Reff 
thus referring to the better constrained inner regions. For the 
fiducial model, we find log gDM(Reff) oc 0.26 log Ly thus ar- 
guing against the universality of this quantity. However, as 
can be seen in Table 4, the slope of the gDAi(Reff)-Ly 
correlation is strongly model dependent with values indicat- 
ing either an increasing behaviour (for NFW models) , a flat 
one (for Burkert + Salpeter) or a decreasing one (for Burk- 
ert + Chabrier). These trends are expected from the analysis 
of column density made above. On the contrary, the scaling 
of g ir {R e ff) does not depend on the halo model because its 
value is estimated from stellar quantities only, while the ef- 
fect of the IMF is simply a systematic rescaling due to the 
higher stellar masses for the Salpeter case. Although the 
slope in Table 4 points towards a significative decrease with 
the luminosity, it is worth stressing that assuming a constant 
value (log g^Reff)} = -9.403 ± 0.01 ± 0.20 (for a Salpeter 
IMF) provides a similarly good fit so that we can not draw 
a definitive answer. 

Any correlation between a DM quantity and a stellar 
one may be the outcome of a hidden interaction between 
the two galactic components. In particular, for the Newto- 
nian accelerations, because of their definition, it is straight- 
forward to show that goM(Reff) = GMdm (Ref /) / R\f / = 
M D M(R eff )/M eff x g*(Ref f ) with M eff = A'L(R eff ) so 
that one can look for a correlation between these two quan- 
tities. For the NFW + Salpeter model, we indeed get : 

log g DM {Ref f ) = 0.21 log g*(Ref f ) - 7.89 . 

with t7i n t = 0.018, the marginalized constraints being : 

a = 0.20j:8:li , P = -7.90±i; 2 g , a int = 0.042t° f . 

Actually, we can recast the above relation in a different way. 
From the definitions of goM(R e ff) and g*(Reff) and the as- 
sumption loggDM(Reff) oc alog g*(R e ff), one easily gets : 

log Mdm (-Re//) oc 2(1 - a) log -R e // + a log M e // 

so that we fit a loglinear relation 

log M DM (R e //) = olm log R eff + (3 M log M eff + ju ■ 



For the best fit relation, we find (1.46, 0.60, 2.79) with ow = 
0.006, while the marginalized constraints (median and 68% 
CL) read 



a M = 1.47 



+0.25 
-0.26 



3 M = 0.61 



+0.22 



JM 



2.54 



+0.41 
-0.26 



(Tint 



= 0.029 



+0.033 



For a = 0.21, we expect <xm — 1-58 in agreement with our 
estimate. On the contrary, j3m is significantly larger than 
a possibly indicating that the ratio gDM{Reff)/g*{R e ff) 
depends on the stellar mass more than expected. However, 
because of the correlation between /3m and (Tint induced by 
the fit, a wrong estimate of the intrinsic scatter may induce 
a bias in the best fit )3m- Since the error bars are quite 
large, determining oi n t is a difficult task as can also be un- 
derstood noting that the best fit Oi nt is formally outside 
the marginalized 68% confidence range (but within the 95% 
one). Nevertheless, the small value of the rms of the resid- 
uals (orms = 0.12) is strong evidence in favour of a very 
tight correlation. We also stress that this scaling relation 
(although with different coefficients) still holds if we change 
the IMF or the halo model. Interestingly, this correlation is 
pretty simila r to th e lu minosity and mass FP d iscussed in 
iBolton et all (|2007l ) and lHvde fc Bernardil (|2009l ). with the 
total MjL ratio found to depend less on the stellar mass 
density than on the effective radius. 

An alternative way to look for the correlation be- 
tween the stellar and DM mass at Re// ma y be 
provided by the DM mass fraction, /dm (Reff) — 

M DM (Ref f) /Wh(R ef f) + M DM (Ref /)]■ Since M DM (R e ;;) 

and M*(i? e //) are correlated, we expect to find a simi- 
lar correlation b etween fpM (R e ff) and mass proxies, such 
as Mi, and Ly (ICappellari et al.ll2006l. IBolton et al.ll2007l . 



iHvde fc Bernardil 120091 . T09. lAuger et al. ||2010hf ). We 
deed find 



log f D M (Reff) = 0.49 log 



log JOM (Reff) = 0.51 log 



Mi, 



10 1 



U 



10 11 u 



-0.86 



0.56 



with (Tint = 0.03 (0.02) and cr rms = 0.13 (0.12) for the first 
(second) case. The marginalized constraints for the /dm - 
Mi, relation are : 



■0A8t° .ll 







-0.86 



+0.08 
0.08 i Pint 



= 0.049 



+0.041 
-0.032 ) 



while for the Jdm - Ly we find : 



(,n+0.09 R „ t -„+0.01 n ma +0.038 

a - 0.51_ .o9 , p = -0.57_o.oi , «»t = 0.028_ . 2i , 

The large error bars on the individual points likely make the 
estimate of (Tint biased, but we nevertheless find clear evi- 
dence for a DM content increasing with both M* and Ly. 
Both these correlations are in very good agreement with 
what is found in T09 for local ETGs for high luminosity 
systems (log Lb ^ 10.4), notwithstanding the different mod- 
els adopted (NFW halo and Sersic stellar profile here vs a 
full isothermal model in T09). However, the slopes we find 
are strongly dependent on the model assumptions. Should 
we use the same NFW model, but the Chabrier IMF, we 
get f DM (Reff) oc Ly 20 and f DM (Reff) oc M°' 16 , which 
are fully consistent with the results we have obtained in 
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ICardone et al. I (|2009T ). where a general galaxy model has 
been fitted to a subsample of SLACS lensea 10 !. 



4 CONCLUSIONS 

Much attention has recently been dedicated to investigating 
whether some correlations can be found among DM quan- 
tities and the stellar ones with contrasting results point- 
ing towards a universal DM column density S (D09, G09) 
or its variation wit h halo mass M200 (B 09) and luminosity 
M* (T09, NRT10. [Auger et al. IfeoiOal L Here we have ad- 
dressed this controversy using a sample of intermediate red- 
shift ETGs using the available data on both the projected 
mass within the Einstein radius and the aperture velocity 
dispersion to fit four different stellar + DM halo models. Mo- 
tivated by recent findings in the literature (|Treu et al.ll2010l . 
NRT10) and the analysis of the virial M/L ratio, we have 
finally chosen a NFW DM halo and a Salpeter IMF which 
we have then used as a reference case for investigating the 
scaling relations of interest. 

Contrary to D09, we find that the column density S, 
evaluated at both the halo characteristic radius r s or the 
stellar effective radius R e ff is not a universal quantity, but 
rather correlates with the luminosity Ly and the stellar and 
halo masses M* and M2oo- Although the slopes of these 
correlations depend on the halo model and IMF, assum- 
ing our reference model, the S(r B )- M200 relation we find 
agrees with the B09 one, with a rather similar slope (0.19 
vs 0.21), but a smaller zeropoint. As a consequence, our 
S(r 3 ) values are smaller than the B09 ones, but also smaller 
than those predicted on the basis of a secondary infall 
model (B 09.lDel Popolo 1 120091') and AC DM N-b ody sim- 
ulations (jMaccio. Dutton & van den Boschl |2008|). We ar- 
gue that this discrepancy is expected considering that these 
studies do not add a stellar component to the galaxy model, 
while here we take this explicitly into account thus decreas- 
ing the DM content. We have also found that the ensemble 
average column density in the central regions is systemat- 
ically larger than the one in spiral galaxies in agreement 
with, e.g., NRT10. This is consistent with mass accretion 
in more massive haloes due to merging of late -type sys- 
tems. As an interesting new result, we have shown that 
a very tight loglinear relation among MnM(R e ff), Reff 
and M e // can be found leading to a DM mass fraction 
which positivel y correlates with both the stellar luminos- 
ity and mass rtCappellari et all I2006J . iBolton et all |200l 
iHvde fe Bernardill2009l . T09V 

The limited mass and luminosity range probed and the 
large errors on the different quantities involved prevent us 
from drawing a definitive answer on the slope and normaliza- 
tion of the above scaling relations. Moreover, a larger dataset 
should also allow us to make a more detailed investigation 
of the impact of halo profiles and IMF assumptions. Should 
these further tests confirm our results, DM scaling relations 
can provide a valuable tool in understanding the physical 
processes which drive galaxy formation and evolution. 
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